Acta Math. Univ. Comenianae 

Vol. (submitted to AMUC), pp. 1- TTul 



1 



SELF-PROPAGATING HIGH TEMPERATURE SYNTHESIS 
(SHS) IN THE HIGH ACTIVATION ENERGY REGIME 

R. MONNEAU and G.S. WEISS 

Abstract. We derive the precise limit of SHS in the high activation energy scal- 
ing suggested by B.J. Matkowksy-G.I. Sivashinsky in 1978 and by A. Bayliss-B.J. 
Matkowksy-A.P. Aldushin in 2002. In the time-increasing case the limit turns out 
to be the Stefan problem for supercooled water with spatially inhomogeneous coef- 
ficients. 

Although the present paper leaves open mathematical questions concerning the 
convergence, our precise form of the limit problem suggest a strikingly simple ex- 
planation for the numerically observed pulsating waves. 



1. Introduction 

The system 

, . d t u - Au = vf(u) 

d t v = -vf(u) , 

where u is the normalized temperature, v is the normalized concentration of the 
reactant and the non-negative nonlinearity / describes the reaction kinetics, is a 
simple but widely used model for solid combustion (i.e. the case of the Lewis num- 
ber being +oo). In particular it is being used to model the industrial process of 
Self-propagating High temperature Synthesis (SHS). In the case of high activation 
energy interesting phenomena like the instability of planar waves, fingering and 
helical waves are observed. 

Since the seventies (and possibly even earlier) it has been argued that the problem 
is for high activation energy related to a Stefan problem describing the freezing of 
supercooled water (see |21| . f 101 p. 57]). In |21| B.J. Matkowsky-G.I. Sivashinsky 
derived a formal singular limit containing a jump condition for the temperature 
on the interface. Later the Stefan problem for supercooled water - the intuitive 
limit - became the basis for numerous papers focusing on stability analysis of Q . 



Received .... 

2000 Mathematics Subject Classification. Primary 80A25, 80A22, 35K55, 35R35. 

Key words and phrases. SHS, solid combustion, singular limit, high activation energy, Stefan 
problem, ill-posed problem, free boundary, hysteresis. 

G.S. Weiss has been partially supported by the Grant-in-Aid 15740100 of the Japanese Min- 
istry of Education and partially supported by a fellowship of the Max Planck Society. Both 
authors thank the Max Planck Institute for Mathematics in the Sciences for the hospitality 
during their stay in Leipzig. 



2 



R. MONNEAU and G.S. WEISS 



fingering, helical waves etc. (see for example [TT | . [12 ] . |9] . [14| . [13 ] . [15 ] . [8] . [T].[2]). 

Surprisingly there are few mathematical results on the subject: In |20| E. Logak- 
V. Loubcau proved existence of a planar wave in one-space dimension and gave a 
rigorous proof for convergence as the activation energy goes to infinity. 
Instability of the planar wave for a special linearization (and high activation en- 
ergy) is due to Qp. 

In the present paper we argue that the SHS system converges to the irreversible 
Stefan problem for supercooled water. As the initial data of the reactant concen- 
tration enters the equation as the activation energy goes to infinity, our result also 
suggests a surprisingly simple explanation for the numerically observed pulsating 
waves (cf. and [2]), namely that they are caused by the spatial inhomogeneity 
ii° (or Y°, respectively) in the below equation and are therefore mathematically 
related to the pulsating waves in [3] 

In the time-increasing case we give a rigorous convergence proof in higher dimen- 
sions. For general initial data in one space-dimension sec our forthcoming paper 

EU 

In the original setting by B.J. Matkowsky-G.I. Sivashinsky |21l equation (2)], 

/ 2 \ dtU N - Au w = (1 - o- N )Ne N v N exp(-N/u N ), 

d t VN = — Ne N Vwexp(— N/un), 

each limit of un > as N — > oo satisfies for (o~N) Ne N CC [0, 1) (for f 
1, N — > oo the limit in this scaling is the solution of the heat equation; cf. Section 
|O an d Theorem P|) 

(3) dtUoo - v°d t x = Aiioo in (0, +oo) x O, 

where v° are the initial data of and 

y(t x) { £ ^' ° SSSUP (o-*) 7 " 00 ^'^ - 1 ' 
\ =1, esssup ( ,t)Woo(-,a;) > 1 > 

and in the time-increasing case, 

)= 0, Uoo(t, x) < 1 , 

e [o,i], Uoo(t,x) = i , 

= 1, Uoo(*, x) > 1 . 

In the SHS system with another scaling and a temperature threshold (see j2] 
p. 109-110]), 
(4) 

dt^N — 

= (1 - a N )NY N exp((iV(l - a N ){6 N - l))/(a N + (1 - o- N )e N ))x { e N> § } , 
d t Y N = -(I - <j n )NY n eM(N(l - <J N )(6 N - l))/(a N + (1 - o- N )e N )) X{ e N> 5} 

where iV(l — ctjv) >> 1, o~n G (0, 1) and 9 G (0, 1), each limit 9oo of On satisfies 
(cf. Section Inl and Theorem WH 

(5) a^oo - Y°d tX = A#oo in (0, +oo) x CI, 
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where Y° are the initial data of and 



e [0, 1], esssup (o^)^ (-,x) < 1 , 
= 1, esssup ( ,t)^oo(-,a;) > 1 , 



and in the time-increasing case, 



r = 0, 0oo(t,») < 1 - 
X (t,x)\ e[0,l], 6 00 (t,x)=l, 

{ = 1, (t,x) > 1 . 



To our knowledge this precise form of the limit problem, i.e. the equation with the 
discontinuous hysteresis term, has not been known. Even in the time-increasing 
case it does not coincide with the formal result in |21| . 

In the case that 9oo (or Uqq, respectively) is increasing in time and v° (or Y°, 
respectively) is constant, our limit problem coincides with the Stefan problem for 
supercooled water, an extensively studied ill-posed problem (for a survey see [3]). 
As it is a forward-backward parabolic equation it is not clear whether one should 
expect uniqueness (see |6j Remark 7.2] for an example of non-uniqueness in a 
related problem). 

On the positive side, much more is known about the Stefan problem for supercooled 
water than the SHS system, e.g. existence of a finger ( |16| 1. instability of the 
finger f |19j ). one-phase solutions those results, when combined with our 

convergence result, suggest that similar properties should be true for the SHS 
system. 

It is interesting to observe that even in the time-increasing case our singular limit 
selects certain solutions of the Stefan problem for supercooled water. For example, 
u(t) = (as — l)X{t<i} + K X{t>i} is f° r each k € (0, 1) a perfectly valid solution of the 
Stefan problem for supercooled water, but, as easily verified, it cannot be obtained 
from the ODE 



Throughout this article R™ will be equipped with the Euclidean inner product x ■ y 
and the induced norm \x\. B r (x) will denote the open n-dimcnsional ball of center 
x , radius r and volume r n w„ . When the center is not specified, it is assumed to 
be 0. 

When considering a set A, \A shall stand for the characteristic function of A, while 
v shall typically denote the outward normal to a given boundary. The operator dt 
will mean the partial derivative of a function in the time direction, A the Laplacian 
in the space variables and C n the n-dimensional Lebesgue measure. 
Finally 1 denotes the parabolic Sobolev space as defined in |18| . 




2. Notation 
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3. Preliminaries 
In what follows, fl is a bounded C -domain in R™ and 

ue& f| w^((o,r)xO) 

Te(0,+oo) 

is a strong solution of the equation 

d t u e (t,x) - Au e (t,x) = -v°(x)d t exp(-i f Q g e (u e {s,x)) ds) , 
1 ' u e (0, •) = u° in fi,V« e -t/ = on (0, +oo) x dfl ; 

here g e is a non-negative function on R satisfying: 

0) <? e is for each e G (0, 1) piecewise continuous with only one possible jump at Zq, 
9e( z o~) = 9ei z o) = in case of a jump, and g e satisfies for each e G (0, 1) and for 
every z E R the bound g e {z) < C £ (l + |z|). 

1) g e /e — » as e — ► on each compact subset of (— oo, 0). 

2) for each compact subset K of (0, +oo) there is ck > such that min(p e , ck) 
ck uniformly on K as e — ► 0. 

The initial data satisfy < u° < C < +oo, u° converges in to w° as e — > 0, 

(u°) c g(o,i) is bounded in L 2 (f2), it is uniformly bounded from below by a constant 
M m ; n , and it converges in _L 1 (f2) to u° as e — > 0. 

Remark 3.1. Assumption 0) guarantees existence of a global strong solution 
for each e G (0, 1). 

4. The High Activation Energy Limit 

Theorem 4.1. XTie family (it e ) ee (o,i) * s / or eac/i T G (0, +oo) precompact in 
L 1 ((0, T) x H), and eac/i limit u of (it e ) eg ( !) as a sequence e rn — » 0, satisfies in 
the sense of distributions the initial-boundary value problem 

(7) d t u - v°d t x = Au in (0, +oo) x fl, 

u(0, ■) = u° + v°H(u°) in f2 , Vu ■ v = on (0, +oo) x 30 , 

wAerc vfi x) I G [ °' 1] ' ° SSSUP (°-'> u( '' x) " ° ' 
[ = 1, esssup , Qt qu(;x) > , 

and H is the maximal monotone graph 

( = 0, 2 < 0, 
£T(zW G[0,1], 2 = 0, 
[ = 1, z > . 

Moreover, \ * s increasing in time and u is a supercaloric function. 

If (u e ) e£ (o l x) satisfies dtu e > in (0,T) x £1, i/ien u is a solution of the Stefan 

problem for supercooled water, i.e. 

dtu — v°dtH(u) = Au in (0, +oo) x . 

Remark 4.2. Note that assumption 1) is only needed to prove the second state- 
ment "If 
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Proof. Step (Uniform Bound from below): 
Since u ( is supcrcaloric, it is bounded from below by the constant u m - ln . 
Step 1 (X 2 ((0,T) x n)-Bound): 

The time-integrated function v e (t, x) := J u e (s, x) ds, satisfies 

(8) d t v £ (t,x) - Av e (t,x) = w t (t,x) +u°(x) 

where w e is a measurable function satisfying < w e < C. Consequently 

f I (d t v e f + \ [ \Vv e \ 2 (T) = [ T [ (w e + u° € )d t v t 
Jo Jn z Jn Jo Jn 

<U T [ (dtv.f + %f (C+\u^ 

implying 



9 / / v ' 9 
z Jo Jn 



(9) / ul <T (C+\u° e \f 

Jo Jn Jn 

Step 2 (L 2 ((0, T) x VL)-Bound for Vmin(it e , M): 
For 

z 2 /2, z < M, 



G M (z) 
and any M 6 N, 



Mz - M 2 /2, z>M, 



T 

G M {ue)~G M {u a e ) + I I |Vmin(u £ ,Af)| 2 
n Jo Jn 

f — min(w £ , M)dt cxp( f g e (u e (s,x)) ds) . 

lo Jn 6 Jo 

As d t exp(— i J Q g t (u e (s, x)) ds) < 0, we know that d t exp(— i J Q <? c (u £ (s, x)) ds) is 

bounded in L°°(0; ^((0, T))), and 

/ -w°min(u e , Af)9 t exp(-i / g e (u e (s,x)) ds) 
o Jn 6 Jo 

< C I sup max(min(u e , M),0) < CM£"(f2) . 
Step 3 (Compactness): 

Let xm : R -* R be a smooth non-increasing function satisfying A(-oo,a/-i) < 
Am < A(-oo,ai) arL< i l°t m be the primitive such that = z for z < M — 1 

and $m < Af. Moreover, let (4>s)se(o,i) be a family of mollificrs. i.e. <j>s G 
C ' (R n ; [0, +oo)) such that J = 1 and supp </>,5 C 5,5(0) . Then, if we extend w £ 
and v® by the value to the whole of (0, +00) x R", we obtain by the homogeneous 
Neumann data of u e that 

d t ($Af(lie) * 4>s) (t,x) 
= (Jxm(u € ) ^AoAu e - v°d t cxp(-j J g e (u e (s,x)) ds)^ * fo^j (t,x) 
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XM{u e )(t, y)[ xn{y)&u e {t,y) 

R" 



1 



v°(y)d t ex P(~- J 3e(« £ ( s : V)) ds )J <Ps(x - y) dy 
4>s(x -y)[ - x'M(u e (t,y))xn{y)\^u e {t,y)\ 2 



R" 



1 



- XM(ut(t,y))v°(y)d t exp(-- I g £ (u e (s,y)) ds) 

e Jo 

+ XM{u e {t, 2/))xn(2/)Vu e (i, J/) • V<fo(x - y) dy . 

Consequently 

/ \dt {§ M {u e ) * 4> s ) I < Ci(fi,C,M,<J,T) 
o JR" 

and 

/ / |V($ M (O*0<5)I < C2(0,M,5,T). 
Jo JR" 

It follows that ($m(w c ) * 05) e g(o i) is for each (M, <5, T) precompact in i 1 ((0, T) x 
R"). 

On the other hand 

f [ |$mW*^-$mW| < C 3 (<5 2 / /|V$ M K)r) 

Jo JR" \ jo jo y 

+ 2(M - Umin )T/:"(B 5 (ar!)) < C!i(C, O, Umi „, M, T) 5 . 

Combining this estimate with the precompactness of (^M(ue)*^s)ee(o,i) we obtain 
that $m(u c ) is for each (M, T) precompact in L 1 ((0, T) x R"). Thus, by a diagonal 
sequence argument, we may take a sequence e rn — > such that <&M(u em ) — > za/ 
a.e. in (0, +oo) x R™ as m — > oo, for every M 6 N. At a.e. point of the set {zm < 
M — 1}, u tm converges to zm- At each point (t, x) of the remainder CIm^nIzm — 
M — 1}, the value u €m (t,x) must for large m (depending on (M, t, x)) be larger 
than M — 2. But that means that on the set C\m£n{ z m > A/ — 1}, the sequence 
( u e m )meN converges a.e. to +oo. It follows that (u Cm ) me N converges a.e. in 
(0, +oo) x fl to a function z : (0, +oo) x £1 — > R U {+oo}. But then, as (ue m ) me N 
is for each T € (0, +oo) bounded in L 2 ((0, T) x f2), (u Cm ) me N converges by Vitali's 
theorem (stating that a.e. convergence and a non-concentration condition in IP 
imply in bounded domains L p -convergence) for each p £ [1, 2) in L p ((0, T) x f2) to 
the weak L 2 -limit u of ("e m ) me N- It follows that £ n+1 (p| MgN {zj\/ > M — 1}) = 
£ n+1 {{u = +oo}) = 0. 

Step 4 (Identification of the Limit Equation in esssup / t \U > 0): 

Let us consider (t, x) G (0, +oo) x SI such that u Cm (s, x) — > u(s, x) for a.e. s <E (0, t) 

and u(-,x) € L 2 ((0,<)). In the case esssup r Qt \u{-,x) > 0, we obtain by Egorov's 

theorem and assumption 2) that exp(— j- J Q g em (u em (s, x)) ds) — > as m — > oo. 
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Step 5 (The case d t u e > 0): 

Let (t,x) be such that u tm (t, x) — > u(t,x) = A < 0: Then by assumption 1), 

/ 1 /"* , / u , x . / max [UmmiA/2] g £ 
cxp( / g €m [u em (s, X)) as) > exp(— t ■ — ) — > 1 as m — > oo . 

Jo ^ m 

Remark 4.3. 7 J for a more general result in one space- dimension see the 
forthcoming paper |22j . 

2) We also obtain a rigorous convergence result in the case of (higher dimensional) 
traveling waves with suitable conditions at infinity. In this case our L (W ' )- 
estimate (Step 2) implies a no- concentration property of the time- derivative. 

5. Applications 

Although the limit equation is an ill-posed problem, the convergence to the limit 
seems to be robust with respect to perturbations of the e-system and the scaling: 
here we mention two examples of different systems leading to the same limit. Other 
examples can be found in mathematical biology (see |17| and ). 

5.1. The Matkowsky-Sivashinsky scaling 

We apply our result to the scaling in [211 equation (2)], i.e. 

. . d t u N - Au N = (1 - a N )Nv N exp(A(l - l/u N )), 

[ ' d t v N = -Nv N exp(A(l - l/u N j), 

where the normalized temperature un and the normalized concentration vm are 
non-negative, (o"at) 7V6 ]\j CC [0, 1) (in the case <jn T 1, N — > oo the limit equation 
in the scaling as it is would be the heat equation, but we could still apply our 
result to itjv/(l — fjv)) and the activation energy N — > oo. 
Setting M m i„ := —1, e := 1/N, u e := — 1 and 

exp((l - l/(z + l))/e),z > -1 
0,z < -1 



and integrating the equation for vn in time, we see that the assumptions of The- 
orem are satisfied and we obtain that each limit itoofoo of un,o~n satisfies 

(11) dtUryo - (1 - (Too)w i9tX = Aiioo in (0, +oo) x fl, 

where 

G [0, 1], csssup ( ,t) u oo{;%) < 1 



{ = 1, esssup ( Oj t)«oo(-,a) > 1 , 
and in the time-increasing case, 

(12) dtUoz - (1 - cr oo )v 9 t i?(w oo ) = Aitoo in (0, +oo) x O, 

Uoc (0, •) = u° + v Q H{u") in n , Viioo • v = on (0, +oo) x 9fi , 

where v° are the initial data of Voq. Moreover, \ is increasing in time and itoo is 
a supercaloric function. 
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5.2. SHS in another scaling with temperature threshold 

Here we consider (cf. p. 109-110]), i.e. 
(13) 

dt^N — A6*at 

= (1 - a N )NY N cxp((iV(l - a N )(0 N - 1))/{(T N + (1 - c7 N )e N )) X{ g N>S} , 
d t Y N = -(1 - a N )NY N exp((2\T(l - a N ){6 N - l))/(a N + (1 - (Tn)O n ))x { o n >§} 

where N(l — <tn) >> l,cr7v G (0,1) and the constant 9 G (0,1) is a threshold 
parameter at which the reaction sets in. 

Setting M min = -1, e := 1/(N(1 - a N )),K(e) := 1 - a N ,u e := 9 N - 1, 



X*) 



exp((z/(K(e)z + l))/e),z > 9-1 
0,z < 0- 1 



and integrating the equation for Y/v in time, we see that the assumptions of The- 
orem ^2 are satisfied and we obtain that each limit Moo of un satisfies 

(14) dtUoo - v°d t x = Auoo in (0, +oo) x Q, 

(t x) { G esssup (o,t) u °o(;x) < 1 , 

' \ =1, esssup ( ,t) u o°{->%) > 1 , 
and in the time-increasing case, 

(15) dtUoo — v°dtH(u oc ) = Auoa in (0, +oo) x Q, 

-"oo(0, ■) = u° + v°H(u°) in n , Vuoo • v = on (0, +oo) x dVL , 

where v° are the initial data of Voq. Moreover, \ is increasing in time and Uoo is 
a supercaloric function. 

6. Open questions 

The most pressing task is of course to study the existence or non-existence of 
"peaking" (cf. Figure of the solution in the negative phase (for the case of one 




t = o 



Figure 1. Is it possible for the solution to have a tiny peak traveling at high speed?. 

space dimension see the forthcoming paper [22]). A related question is whether 
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(w e ) ee (o,i) is bounded in L°° in the case of uniformly bounded initial data. Al- 
though this seems obvious, it is not obvious how to prevent concentration close to 
the interface. 

Uniqueness for the limit problem (the irreversible Stefan problem for supercooled 
water) in general seems unlikely. One might however ask whether time-global 
uniqueness holds in the case that u is strictly increasing in the xi-dircction. By 
the result in [7j for the ill-posed Hclc-Shaw problem, time-local uniqueness is likely 
to be true here, too. 

Acknowledgment. We thank Stcphan Luckhaus. Mayan Mimura, Stefan 
Miillcr, and Juan J.L. Velazquez for discussions. 
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